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摘 要 : 快速 获取 区 域 土壤 盐 涡 化 程度 信息 ,对 


盐 


渍 化 治理 与 生态 环境 保护 具有 重要 意义 。 以 银 


川 平 原 为 研究 区 ,以 盐分 影响 因子 和 盐分 指数 分 别 作为 输入 参数 ,建立 支持 向 量 机 (SVM) ,BP 神 经 
网 络 (BPNN ) 和 贝 叶 斯 神经 网 络 (BNN)3 种 土壤 盐分 预测 模型 ,选取 最 佳 模型 进行 研究 区 不 同 深度 
的 土壤 盐 渍 化 预测 。 结 果 表 明 :(1) 0~20 cm 土壤 盐分 预测 模型 中 基于 影响 因子 变量 组 的 BNN 模 
型 效果 最 佳 RE ABR’) HW 0.618, 247 HIRE (RMSE) X 2.986; 20~40 cm 土壤 盐分 预测 模型 中 基 
于 盐分 指数 变量 组 的 BNN 模型 效果 最 佳 , 尼 为 0.651,RMSE 为 1.947; 综合 对 比 下 ,BNN 模型 的 预测 
效果 最 好 ,可 用 于 研究 区 土壤 盐 渍 化 预测 。(2) 银川 平原 主要 是 以 非 盐 渍 化 和 轻 度 盐 渍 化 为 主 ,0~ 


20 cm 土壤 重度 盐 渍 化 及 盐 土 共 占 总 


积 的 11.59% ,20~40 cm 土壤 重度 盐 渍 化 及 盐 土 共 占 总 面积 


的 7.04% ,20~40 cm 土壤 盐 渍 化 程度 较 0~20 cm 土壤 盐 渍 化 轻 。 
关键 词 : 机 器 学 习 ; 土壤 盐分 预测 ; 贝 叶 斯 神经 网 络 ; 银川 平原 
文章 编号 : 1000 - 6060(2023)01 - 0103 - 12(0103 ~ 0114) 


土壤 是 人 类 进行 社会 生产 生活 的 基本 条 件 , 土 
壤 退 化 牵制 着 农业 的 发 展 ,监测 土壤 性 质 是 一 个 重 
要 过 程 ,对 如 何 可 持续 利用 起 着 重要 作用 。 全 球 


—= 


进行 大 范围 区 域 的 盐分 监测 且 优 于 多 元 线性 回归 
模型 和 地 理 加 权 回 归 模 型 。 杨 练兵 等 "在 BPNN 模 


型 基础 上 使 用 遗传 算法 优化 输入 参数 子 集 和 隐 含 


六 大 洲 的 土壤 中 均 有 盐 渍 化 现象 发 生 ( 除 南极 洲 
未 有 调查 数据 外 ) ,并 以 每 年 约 1.5xl0' km 的 速度 
增加 玉 。 盐 分 预测 指 的 是 将 某 些 指标 与 其 相对 应 的 
地 表 信 息 间 构建 关系 模型 ,众多 学 者 对 盐 渍 化 预测 
展开 了 相关 人 研究 ~ 。 徐 红 涛 等 " 姜 红 等 "对 新 疆 
不 同 区 域 进 行 土壤 盐 渍 化 预测 模型 的 构建 ,研究 表 
明 机 器 学 习 的 方法 使 得 模型 的 预测 精度 高 于 传统 
的 回归 法 预测 模型 。 马 国 林 等 "基于 Sentinel-2A 影 
像 , 使 用 XGBoost 选 择 有 效 变量 后 结合 机 器 学 习 进 
行 土壤 盐分 反 演 模型 训练 ,模型 预测 效果 明显 提 
高 ,使 得 盐 渍 化 制图 精度 得 到 提升 。 杨 厚 翔 等 ” 选 
取 多 个 影响 盐 碱 化 的 因子 ,借助 BP 网络 模型 对 黑龙 
江 省 地 区 进行 盐 碱 化 危险 度 监测 ,整体 准确 度 高 于 
95%。 刘 全 明 等 利用 微波 雷达 数据 对 盐分 进行 反 
演 , 所 得 BP 神经 网 络 (BPNN ) 模 型 能 够 在 短 时 间 内 
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层 神经 元 数量 并 对 模型 的 初始 权重 进行 优化 ,确定 
了 适用 于 区 域 的 反 演 模型 ,使 得 反 演 结果 中 土壤 盐 
分 含量 的 均 质 性 最 好 。BPNN 模型 的 使 用 ,提升 了 
人 们 对 于 土壤 盐 渍 化 变化 及 预测 的 精度 ,丰富 了 盐 
渍 化 研究 在 时 空 方面 的 意义 。 

人 工 神经 网 络 在 预测 方面 已 经 有 了 较 多 的 应 
用 成 果 , 章 龙 管 等 将 贝 叶 斯 网 络 应 用 于 施工 风险 
管理 中 ,结果 证 明 该 方法 对 于 施工 风险 及 发 展 趋势 
能 够 有 效 预测 ,可 为 现场 施工 防范 管理 提供 理论 指 
导 。 毕 春光 等 所 对 玉米 病害 进行 预警 研究 ,证明 贝 
叶 斯 神经 网 络 (BNN ) 模 型 较 BP 神 经 模型 的 精准 度 
提升 了 5.49%。 将 BNN 模 型 应 用 到 土壤 盐 渍 化 预测 
中 ,探讨 不 同 模型 对 土壤 盐分 的 拟 合 效果 ,本文 以 
银川 平原 为 研究 区 ,基于 遥感 指数 和 盐 渍 化 影响 因 
子 ,结合 实测 土壤 盐分 数据 构建 预测 模型 , 选 出 最 
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优 模型 ,预测 研究 区 不 同 土壤 深度 的 盐分 分 布 ,为 
银川 平原 土壤 监测 和 盐 渍 化 防治 提供 理论 支持 。 


1 研究 区 概况 


银川 平原 (37.83°~39.38°N, 104.28°~107.65°E) 
西北 、 东 三 面 被 腾 格 里 沙漠 .马兰 布 和 沙漠 ES 
素 沙 漠 环 绕 , 南 部 与 黄土 高 原 接壤 。 银 川 平原 属于 
温带 干旱 区 ,多 年 平均 气温 5~9 %C ,年 平均 降水 量 约 
200 mm, Z$ K Lt 24 10:1, 42 H AR AY He 2500~3100 h, 
年 太阳 辐射 5800x105~6100x105J.m2。 银 川 平 原 主 
要 是 由 冲 洪 积 平原 构成 ,地 势 西南 高 东北 低 。 黄 河 
贯穿 银川 平原 ,优越 的 地 理 位 置 使 得 平原 能 够 充分 
发 展 水 利 灌溉 ,灌溉 渠 系 众多 , 农 牧 业 发 达 '"" ,主要 
作物 有 水 稻 小麦、 玉米 。 灌 洲 土 和 草 甸 土 是 主要 
的 农业 土壤 ,熟化 程度 较 高 。 地 下 水 矿 化 度 为 0~6 
g L ', 平 原 地 区 地 下 水 埋 深 较 浅 ,土壤 盐分 累积 


2 数据 与 方法 


2.1 土壤 样本 采集 与 盐分 测定 

本 文大 致 以 5 kmx5 km 的 格 网 设置 采样 点 (图 
1), F 2021 463 A 7—28 日 采样 。 采 样 时 每 个 样 点 
采用 梅花 形 采样 法 ,在 采样 点 30 mx30 m 范 围 内 采 
集 0~20 cm 和 20~40 cm 土壤 混 匀 ,运用 四 分 法 留 约 
500 g 土 壤 装 进 密封 袋 带 回 ,并 用 GPS 记录 采样 点 的 
坐标 信息 及 周围 环境 信息 , 共 采 集 197 个 样本 点 。 
采集 的 土壤 样本 剔除 杂质 后 自然 风干 ,研磨 后 过 2 
mm 筛子 ,以 水 土 比 S$:1 配 置 提取 液 二 ,每 个 样品 配 
置 3 组 作为 对 照 组 ,每 个 样 点 取 3 组 数值 的 平均 值 
作为 样 点 最 终 pH 和 电导 率 , 运 用 电导 率 法 "计算 
土壤 样本 的 含 盐 量 。 

S=(EC+41.2543)x 5 + 2120.76 (1) 

式 中 :3 为 土壤 含 盐 量 (g'kg );EC 为 测量 所 得 土壤 
提取 液 电导 率 (S.m-”)。 

对 测量 结果 进行 分 析 ,删除 异常 样本 点 后 获得 
166 个 样 点 供用 。 将 样本 的 盐 渍 化 程度 分 为 $ 级 一 ， 
样本 的 统计 性 描述 如 表 1 所 示 。0~20 cm 土壤 样本 
中 共有 91 个 盐 渍 化 样本 , 约 占 总 样本 的 54.829% 。 
20~40 cm 土壤 样本 共有 65 个 盐 渍 化 样本 , 约 占 总 样 
本 的 39.16%。 各 层 样本 土壤 含 盐 量 总 体 变异 系数 
均 超 过 132% ,变异 性 较 强 ,表明 样本 离散 程度 较 


图 1 银川 平原 采样 点 分 布 


Fig. 1 Distribution of sampling points in Yinchuan Plain 


高 ,具有 普 适 性 。 

将 土壤 样本 按 含 盐 量 由 低 到 高 排序 ,依据 模型 
训练 集 和 结果 验证 集 2:1 的 比例 等 间隔 选择 110 个 
样本 进行 模型 训练 ,56 个 样本 作为 验证 集 用 于 模型 
ue, 

2.2 遥感 影像 获取 与 处 理 

研究 选用 的 Landsat 8 OLI 遥 感 影像 来 源 于 美国 
地 质 调查 局 官网 (https://earthexplorer.usgs.gov/) , 行 
列 号 为 129/033 和 129/034, 空 间 分 辨 率 为 30 m, 成 
像 时 间 为 2021 年 3 月 17 日 ,与 采样 时 间 同 期 。 利 用 
ENVI5.3 对 所 选 的 影像 进行 大 气 校正 等 预 处 理 流 程 
后 ,进行 各 类 盐分 指数 的 计算 和 提取 ” ,如 表 2 
所 示 。 

2.3 土壤 盐 渍 化 预测 模型 

2.3.1 盐 渍 化 影响 因子 组 参数 选择 ”依据 已 有 学 者 
对 土壤 盐 渍 化 成 因 条 件 的 相关 人 研究”, 结合 银川 
平原 实际 情况 ,选择 8 个 环境 变量 和 3 个 人 为 活动 
变量 ( 表 3)。 其 中 ,土壤 pH 和 含水 率 结果 由 采样 点 
数据 计算 后 进行 Kriging 空间 插值 而 得 ,设置 输出 的 
分 辩 率 为 30 m; 地 下 水 数据 来 源 参 考 文献 [17] ,将 
数据 处 理 为 30 mm 分 辨 率 的 栅 格 数据 ;银川 平原 数字 
高 程 模型 从 地 理 空间 数据 云 平台 (http://www.gs- 
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表 1 银川 平原 土壤 样本 统计 
Tab. 1 Statistics of soil sample in Yinchuan Plain 
深度 /em 程度 Ab lible ke! 样本 数 ares aaa a as ba a 
0~20 非 盐 涡 化 <1 75 0.510 0.983 0.097 41.295 172.327 
轻 度 盐 渍 化 1~2 32 1.522 1.946 1.058 18.823 
中 度 盐 渍 化 2~4 20 2.819 3.667 2.003 20.119 
重度 盐 渍 化 4~6 12 4.734 5.803 4.122 11.538 
tht >6 27 15.104 43.725 6.676 59.548 
20~40 非 盐 涡 化 <1 101 0.566 0.997 0.308 30.428 132.278 
轻 度 盐 渍 化 1~2 34 1.471 1.979 1.020 17.701 
中 度 盐 渍 化 2~4 23 2.674 3.752 2.015 19.605 
重度 盐 渍 化 4~6 4 4.387 5.276 4.038 13.545 
tht >6 4 9.872 17.929 6.274 55.894 


表 2 盐分 指数 计算 公式 


Tab. 2 Calculation formulas of salt indices 


盐分 指数 计算 公式 参考 文献 
SI-T 100(Red-NIR) SURO 
NDSI Red-NIR)/(Red+NIR) Amal 22 
Sl Blue/Red Sahana 等 1 
S2 Blue—Red)/(Blue+Red) Nguyen 等 中 
S3 GreenxRed)/Blue Nguyen £% 
S5 BluexRed)/Green PIN Ag 
S6 RedxNIR)/Green 孙 亚 楠 等 中 
SIl Green x Red HE pe Fa ee 
S13 Green’ + Red? 赵 巧 珍 等 中 


注 :Blue ,Green Red NIR 分 别 为 蓝 波段 (b2) . 绿 波 段 (b3) ` 红 波段 
(b4)、 近 红外 波段 (b5 ) 的 光谱 反射 率 。 下 同 。 


cloud.cn/) 获 取 , 空 间 分 辩 率 为 30 m; 地 表 温 度 .增强 
植被 指数 和 水 体 指数 通过 Landsat 8 OLI 影 像 在 EN- 
VI 5.3 软件 中 计算 获得 ;土地 利用 数据 从 中 国 科学 
院 资 源 环境 科学 与 数据 中 心 (http://www.resdc.cn/) 
获取 ,强度 指数 采用 庄 大 方 ” 利 用 分 级 方法 ,土地 
利用 强度 依据 不 同 的 土地 类 型 分 别 确定 为 不 同 的 
值 ,未 利用 地 、 草 地 、 林 地 、 耕 地、 居民 用 地 水体 分 
别 赋 值 为 1.2、2、3、4、1。 

2.3.2 数据 处 理 为 消除 数据 获取 过 程 中 人 为 或 者 
自然 因素 导致 数据 不 标准 的 问题 , 先 将 数据 进行 归 
一 化 处 理 以 消除 不 同 指标 之 间 维 度 的 影响 。 归 一 
化 公式 为 : 


A-A,, 
A= a (1) 


max min 


NP: 4 为 样本 数据 的 最 大 值 ; 4,, 为 样本 数据 
的 最 小 值 ; 4 为 预警 因子 数据 集 ; A, 为 归 一 化 后 的 
数据 。 

2.3.3 预测 模型 构建 

(1) 支持 向 量 机 (SVM ) 

SVM 模型 基于 结构 风险 最 小 原理 ,依靠 有 限 的 
样本 来 检索 全 局 最 优 解 ,拥有 对 未 知 点 有 较 好 的 泛 
化 效果 的 优点 呈 ,规避 小 样本 学 习 和 局 部 极 值 出 现 
的 缺点 。 本 文选 用 RBF(Radial basis function) 为 
SVM 中 核 函 数 类 型 ,惩罚 参数 (ec) 与 核 参量 (g) 由 样 
本 进行 测试 计算 得 到 。 

(2) BP 神经 网 络 (BPNN) 

BPNN 模 型 基于 反 疝 传播 误差 的 方式 对 数据 集 
进行 训练 ,达到 误差 最 小 化 的 目的 ”"。 该 模型 拥有 
较 强 的 自 适应 和 自学 习 能 力 , 可 精确 通 近 任意 的 非 
线性 关系 ?101。 本 文 基于 MATLAB 2020 软件 ,将 盐分 
指数 和 影响 因子 作为 模型 的 输入 数据 ,实测 盐分 数 
据 作为 输出 数据 ,创建 研究 区 土壤 盐分 的 BPNN il 
测 模型 。 经 多 次 训练 ,在 参数 一 致 的 情况 下 ,将 模 
型 的 输出 值 和 实际 值 进行 拟 合 , 以 误差 作为 选择 隐 
含 层 节点 的 标准 ,确定 最 佳 训练 效果 的 隐 含 层 节 
点 。 模 型 训练 目标 的 最 小 误差 设置 为 0.00001 ,训练 
次 数 设置 为 1000 次 。 

(3) 贝 叶 斯 神经 网 络 (BNN) 

BNN 模 型 通过 为 神经 网 络 的 权重 分 配 引 入 随 
机 性 进行 正则 化 ,也 相当 于 任意 权重 的 多 种 神经 网 
络 的 预测 ”31。 与 传统 BPNN 模 型 不 同 的 是 贝 叶 斯 
估计 得 到 的 是 节点 权重 的 后 验 分 布 概率 P(wlD) ,并 
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表 3 土壤 盐 渍 化 影响 因子 


Tab.3 Influence factors of soil salinization 


变量 类 型 参数 名 称 描述 
环境 变量 土壤 pH 土壤 酸 碱 性 
土壤 含水 率 /% 土壤 绝对 含水 量 
高 程 /m 地 区 的 海拔 高 度 
增强 植被 指数 农作物 长 势 和 营养 信息 
水 体 指数 地 区 水 体 信息 
地 表 温度 /*C 也 区 土壤 温度 
地 下 水 埋 深 /m 也 区 洪水 的 埋藏 深度 
地 下 水 矿 化 度 /g:L" 区 域 潜水 中 可 溶性 盐 类 的 质量 
人 为 活动 变量 土地 利用 强度 土地 资源 的 利用 效率 
人 口 密度 /人 :km” 每 平方 千 米 的 人 口 平 均 数 
生产 总 值 /元 生产 活动 的 最 终 成 果 


不 是 固定 的 值 ,因此 可 以 为 神经 网 络 的 预测 添加 随 
机 性 。 由 给 定数 据 输 出 预测 值 的 分 布 ,公式 为 : 

P(nlm) = E pap |P(nlm,w)| (2) 
式 中 :给 定数 据 、 预 测 值 分 别 为 m 和 nn;w 为 每 个 节点 
的 权重 ;D 为 训练 集 。 下 同 。 

所 得 到 的 预测 值 y 是 P(wID) 中 的 所 有 可 能 的 预 
测 值 的 期 望 。 依 据 贝 叶 斯 理论 ,样本 数据 决定 P(D)。 
P(w,D) £ P(Dw)P(w) (3) 

P(D) P(D) 
根据 参考 文献 [34] ,获得 一 个 高 斯 完 验 的 Pw) : 
P(w)=TIrN(w0,0?)+ 0 -nm)N(w0,0o3) (4) 


式 中 : w 为 网 络 的 第 j 个 权重 ; oil 、o 为 先 验 的 第 
一 个 和 第 二 个 混合 成 分 的 方差 , 且 o>; 
N(wi0, 中 为 wj 以 均值 为 0 评价 的 高 斯 密度 。 
依据 所 给 的 先 验 分 布 ,运用 贝 叶 斯 定理 得 到 后 
验 概率 ,若菜 参数 能 够 使 后 验 概率 达到 最 大 化 , 则 
该 参数 为 模型 的 最 佳 参 数 。BNN 模型 能 够 克服 过 
拟 合 及 出 现 局 部 最 小 化 的 现象 。 
2.3.4 模型 评价 指标 为 量化 土壤 盐分 预测 模型 的 
预测 效果 ,本 文选 择 决定 系数 (RR) 和 均 方 根 误差 
(CRMSE)2 个 常用 指标 来 进行 综合 评价 。 尼 的 值 越 
高 ,RMSE 的 值 越 低 , 说 明 模型 拟 合 效 果 越 好 。 


3 结果 与 分 析 


3.1 银川 平原 土壤 盐 渍 化 空间 分 布 特征 分 析 
[EH ArcGIS 的 空间 分 析 工 具 对 166 个 采样 点 的 


P(wlD) = 


土壤 含 盐 量 进行 插值 ,根据 盐 渍 化 程度 划分 为 5 个 
等 级 ,得 到 银川 平原 不 同 深度 土壤 含 盐 量 的 空间 分 
布 特征 (图 2)。 由 图 2 可 见 ,0~20 cm 重度 盐 渍 化 土 
壤 及 盐 土 主要 分 布 在 银川 平原 北部 地 区 ,南部 地 区 
盐 涡 化 程度 较 轻 。20~40 cm 主要 以 非 盐 渍 化 和 轻 
度 盐 渍 化 土壤 为 主 ,银川 平原 西北 部 土壤 盐 渍 化 现 
象 严重 ,部 分 区 域 土壤 含 盐 量 大 于 6 eke. BARE 
研究 区 0~20 cm 和 20~40 cm 土壤 盐 渍 化 呈现 出 南 低 
北 高 的 趋势 ,重度 盐 渍 化 土壤 和 盐 土 主要 分 布 在 石 
中 山地 区 。 

3.2 不 同 预测 模型 对 比分 析 

3.2.1 相关 性 分 析 对 盐 渍 化 影响 因子 .盐分 指数 
与 实测 土壤 盐分 进行 相关 性 分 析 ,所 得 结果 如 表 4 
所 示 。 由 表 4 可 见 ,影响 因子 中 ,0~20 cm 土壤 除 土 
壤 pH .土壤 含水 率 、 人 口 密度 外 ,其 余 因 子 与 含 盐 量 
具有 较 强 的 显著 性 ;20~40 cm 土壤 含 盐 量 与 土地 利 
用 强度 .增强 植被 指数 ,水 体 指 数 `. 地 表 温 度 .地 下 
水 埋 深 、 地 下 水 矿 化 度 在 0.01 水 平 上 具有 极 强 的 显 
著 性 。 盐 分 指数 中 ,0~20 cm 土壤 售 盐 量 与 S1、52、 
S3、S6、SI2、NDSI、SI-T 的 相关 性 均 通过 0.05 检验 。 
20~40 cm 土壤 含 盐 量 与 S1.S2 .S4 .S5 、SI3、NDSI、 
SIT 的 相关 性 均 通 过 0.05 检 验 。 故 0~20 cm 土壤 
择 通 过 0.01 显著 性 水 平 检验 的 影响 因子 和 盐分 指 
数 分 别 作 为 输入 层 参数 进行 建 模 ,20~40 cm 土壤 选 
择 通 过 0.05 显著 性 水 平 检验 的 影响 因子 和 盐分 指 
数 分 别 作为 输入 层 参数 进行 建 模 。 

3.2.2 模型 构建 与 对 比分 析 选择 SVM、BPNN 和 
BNN 模型 ,以 土壤 含 盐 量 为 输出 层 ,相关 性 强 的 影 
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(b) 20~40 cm 


图 例 
土壤 盐 渍 化 程度 
mm JERE 
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m it 
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2 土壤 盐分 空间 分 布 
Fig.2 Spatial distribution of soil salinity 


响 因 子 和 盐分 指数 为 输入 层 数 据 ,对 0~20 cm Fil 20~ 
40 cm 土壤 盐分 含量 进行 模拟 ,结果 如 表 5 所 示 。 以 
验证 集 样本 对 模型 进行 验证 ,结果 如 图 3。 由 表 5 可 
见 ,0~20 cm 土壤 影响 因子 变量 组 的 3 个 模型 尼 均 大 
于 0.76, 其 中 BNN 模 型 建 模 效果 最 好 ,BPNN 模 型 次 
之 ,SVM 模 型 最 差 。 由 验证 结果 表明 ,BNN 模 型 的 
验证 效果 与 SVM 模型 和 BPNN 模型 相 比 结果 最 佳 ， 
其 模拟 值 与 实测 值 的 尼 值 最 大 (图 3)。 盐 分 指数 组 
的 尺 整体 较 低 ,但 均 通 过 0.01 显著 性 检验 。 其 中 ， 
该 分 组 中 BNN 模 型 的 验证 集 RSME 为 最 低 ,BPNN 
模型 拟 合 误差 最 大 。 对 比 6 个 模型 的 验证 效果 ,发 
现 基 于 影响 因子 组 建立 的 BNN 模 型 是 0~20 cm 土壤 
盐分 预测 最 佳 模型 , 故 选 择 此 模型 对 银川 平原 0~20 
om 土壤 盐分 进行 预测 。 

20~40 cm 土壤 盐分 预测 模型 中 ( 表 5) ,影响 
子 组 模型 验证 集 的 RMSE 均 大 于 2,BNN 模型 的 忆 
大 于 SVM 和 BPNN 模型 ,BNN 模型 的 RMSE 小 于 
SVM 和 了 BPNN 模型 ,因此 BNN 模型 检验 结果 效果 优 
FSVM 和 BPNN 模 型 (图 3)。 盐 分 指数 组 验证 集 的 
尼 相 比 于 影响 因子 组 有 明显 提升 ,盐分 指数 组 BNN 


模型 是 3 个 模型 中 误差 最 小 、 拟 合 程度 最 佳 的 模型 ， 
可 以 将 该 模型 用 于 研究 区 20~40 cm 土壤 盐分 预测 。 
3.3 银川 平原 土壤 盐分 预测 

利用 经 过 训练 的 BNN 模型 对 整个 银川 平原 0~ 
20 cm 、20~40 cm 土壤 进行 盐分 预测 ,而 后 按照 盐 渍 
化 程度 对 结果 划分 等 级 ,得 到 银川 平原 土壤 盐分 预 
测 等 级 图 (图 4) ,由 图 4 可 知 ,0~20 cm 土壤 中 盐 土 
主要 分 布 在 银川 平原 北部 的 大 武 口 区 和 平 罗 县 , 重 
度 盐 渍 化 和 中 度 盐 渍 化 土壤 主要 分 布 在 盐 土 四 周 
及 平原 东部 地 区 , 非 盐 渍 化 和 轻 度 盐 渍 化 土壤 主要 
分 布 在 平原 南部 地 区 。0~20 cm 土壤 主要 是 以 非 盐 
涡 化 土壤 为 主 , 盐 圭 分布 于 大 武 口 区 及 平 罗 县 西部 
地 区 ,重度 盐 涡 化 土壤 零星 分 布 在 银川 平原 内 。 根 
据 像 元 计算 得 到 ,0~20 em 土壤 主要 是 以 非 盐 渍 化 
土壤 为 主 , 占 银川 平原 总 面积 的 33.29% ,其 次 为 轻 
度 盐 渍 化 土壤 ,面积 占 比 为 31.86% ,中 度 盐 渍 化 和 
重度 盐 渍 化 土壤 面积 占 比 分 别 为 23.26% 和 8.07% ， 
盐 土 面积 占 比 (3.52% ) AE 5%. 20~40 cm 土壤 中 
非 盐 涡 化 土壤 和 轻 度 盐 渍 化 土壤 面积 占 比 分 别 为 
50.54% 和 30.07% , 中 度 盐 渍 化 土壤 面积 占 比 为 
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表 4 模型 参数 与 土壤 盐分 的 相关 性 


Tab. 4 Correlation between model parameters and soil salinity 


参数 类 型 参数 与 0~20 cm 土壤 相关 系数 (7) 与 20~40 cm 土壤 相关 系数 (7,) 

盐 渍 化 影响 因子 土壤 pH 0.034 0.030 
土壤 含水 率 0.103 0.102 
高 程 -0.253” -0.114 
土地 利用 强度 -0.403" -0.255" 
增强 植被 指数 -0.580” -0.311” 
水 体 指数 0.511” 0.271" 
地 表 温 度 -0.445" -0.119" 
地 下 水 埋 深 -0.477" -0.204” 
地 下 水 矿 化 度 0.768” 0.254” 
人 口 密度 -0.023 0.048 
生产 总 值 -0.173° -0.025 

盐分 指数 S1 0.492™ 0.402" 
S2 0.480" 0.388” 
S3 -0.230" -0.054 
S4 0.049 0.210” 
S5 0.067 0.229" 
s6 -0.309” -0.093 
SII -0.032 0.136 
SI2 -0.130 0.069 
SI3 -0.043 0.127" 
NDSI 0.325" 0.204” 
SLT 0.347" 0.182" 


注 :*、** 分 别 表示 P<0.05 .P<0.01 水 平 上 的 显著 性 检验 。 


表 5 基于 不 同 变量 组 的 机 器 学 习 模 型 


Tab.S Machine learning models based on different variable groups 


nem 模型 类 别 建 模 集 (n=110) 验证 集 (n=56) 

R RMSE R RSME 

影响 因子 组 SVM 0.766 3.397 0.542 3.266 
(0~20 cm) BPNN 0.773 2.638 0.608 3.081 
BNN 0.797 2.751 0.618 2.986 

盐分 指数 组 SVM 0.319 5.945 0.213 4.155 
(0~20 em) BPNN 0.311 5.878 0.206 4.405 
BN 0.336 5.715 0.250 3.890 

影响 因子 组 SVM 0.368 1.327 0.221 2.401 
(20~40 em) BPNN 0.314 1.353 0.331 2.221 
BN 0.425 1.076 0.486 2.026 

盐分 指数 组 SVM 0.440 1.032 0.436 2.079 
(20~40 cm) BPNN 0.348 1.086 0.545 2.033 
BNN 0.442 1.006 0.651 1.947 


注 :SVM BPNN BNN 分 别 为 支持 向 量 机 、BP 神 经 网 络 和 贝 叶 斯 神经 网 络 ; 尼 为 决定 系数 ;RMSE 为 均 方 根 误差 。 


12.35% ,重度 盐 清 化 土壤 和 盐 土 面积 占 比 (分 别 为 人 性 ,将 预测 结果 与 实测 盐分 插值 结果 进行 栅 格 对 
4.89% .2.15% ) 均 小 于 5$%。 LE ,统计 相同 等 级 像 元 数据 ( 表 6)。 由 表 6 可 知 ,0~ 
为 进一步 验证 预测 结果 与 实验 室 结果 的 一 致 20 cm 和 20~40 cm 土壤 盐 渍 化 预测 所 得 的 盐 渍 化 程 
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(c) 0~20 cm 影响 
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图 3 土壤 盐分 与 模型 预测 值 的 散 点 图 


Fig. 3 Scatter plot of soil salinity and model predictions value 


度 与 实测 含 盐 量 插值 相同 等 级 像 元 个 数 的 相同 率 
均 大 于 73% ,表明 预测 模型 的 结果 与 实测 数据 间 具 
有 较 好 的 一 致 性 ,所 选 模型 可 达到 检验 要 求 并 用 于 
研究 区 的 土壤 盐 渍 化 预测 。 


4 讨论 


使 用 3 种 机 器 学 习 的 方法 建立 了 对 于 不 同 深度 
的 土壤 盐分 预测 模型 ,发 现 不 同 深度 的 土壤 对 同 种 
变量 在 相同 模型 下 的 建 模 效果 差异 较 大 ,对 所 有 结 


果 进 行 精 度 评价 后 ,选择 适用 不 同 深度 土壤 的 最 佳 
模型 进行 土壤 盐分 预测 。 

本 文 对 多 个 影响 因子 与 实测 值 进行 Pearson 相 
关 性 分 析 ,发 现 敏感 因子 为 地 下 水 矿 化 度 .地 下 水 
埋 深 .土地 利用 强度 等 多 个 因子 ,与 杨 思 存 等 ”、 苏 
春 利 等 “所 得 结论 相符 ,后 续 可 用 于 挖掘 因素 与 土 
壤 盐 分 的 隐 含 信息 。0~20 cm 土壤 盐分 与 各 影响 因 
子 的 相关 性 明显 高 于 20~40 cm 土壤 ,这 与 刘 继 龙 
等 ”研究 结果 一 致 ,可 能 是 土壤 盐分 受到 外 界 影 响 
程度 随 深 度 不 同 而 改变 ,人 们 所 进行 的 农业 活动 主 


202302.00266v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


110 FFE WU 46 卷 
(a) 0~20 cm (b) 20~40 cm 
N N 


人 


图 例 
= 水 体 

mm 城镇 用 地 
土壤 盐 涡 化 程度 
me SEED tL 
me 轻 度 盐 渍 化 
国 | 中 度 盐 渍 化 
国 | 重度 盐 渍 化 
me fet 

0 10km 


—— 


图 4 土壤 合 盐 量 
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Fig. 4 Grade distribution of soil salinity inversion 


表 6 实测 插值 与 预测 像 元 统计 


Tab.6 Measured interpolation and predictive pixel statistics 


深度 /cm 盐 渍 化 程度 实测 插值 像 元 数 预测 和 实测 插值 相同 等 级 个 数 相同 率 /% 

0~20 非 盐 渍 化 3922460 3221641 82.13 
轻 度 盐 渍 化 1925470 1409137 73.18 
中 度 盐 渍 化 779011 648843 83.29 
重度 盐 渍 化 764161 602002 78.78 
HE 507095 434657 85.72 

20~40 非 盐 渍 化 4279527 3613169 84.43 
轻 度 盐 渍 化 2413166 2046062 84.79 
中 度 盐 渍 化 722601 538546 74.53 
重度 盐 渍 化 311423 274697 88.21 
盐 土 171493 144933 84.51 


要 是 在 0~20 cm 土壤 ,对 20~40 cm 土壤 的 干预 性 
较 小 。 

边 慧 芹 等 ”\ 何 宝 忠 等 ”对 于 不 同 地 区 的 盐 汗 
化 进行 研究 ,表明 盐分 指数 在 监测 应 用 方面 具有 一 
定 的 优势 ,选择 多 种 盐分 指数 与 实测 盐分 进行 相关 
性 分 析 ,选取 相关 性 较 高 的 盐分 指数 作为 模型 参数 
进行 预测 模型 训练 。 本 文采 样 时 正 值 春季 ,气候 干 
IR ,蒸发 强烈 ,土壤 盐分 随 水 分 的 蒸发 而 在 地 表 集 


聚 ,使 0~20 cm 土壤 含 盐 量 高 于 20~40 om 土壤 ,造成 
盐分 指数 与 两 层 土壤 盐分 的 相关 性 差异 。 土 壤 是 
一 个 复杂 的 综合 体 ,常规 的 线性 回归 模型 难以 达到 
较 好 的 精度 ,机 器 学 习 能 够 较 好 地 克服 这 个 难题 ， 
许多 研究 者 将 机 器 学 习 用 于 预测 盐分 ,所 得 模型 精 
度 高 于 线性 模型 ss BNN 算法 可 运用 较 少 的 数据 
得 到 稳定 的 模型 和 各 层 参 数 的 分 布 ,达到 较 好 的 泛 
化 目的 ,可 为 盐分 预测 提供 新 思路 。 本 文 基于 2 个 
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变量 组 ,对 3 种 不 同 机 器 学 习 方法 进行 建 模 比 较 , 发 
现 BNN 模型 精度 优 于 BPNN 模型 和 SVM 模型 。 且 
0~20 cm 土壤 盐分 预测 整体 效果 优 于 20~40 cm 深度 
的 模型 效果 ,与 杨 宁 等 ”所 得 结论 基本 一 致 。 

本 文通 过 对 比 不 同 深度 和 不 同 输入 参数 下 多 
个 模型 的 预测 效果 ,选择 最 佳 模 型 进行 预测 。 同 一 
地 区 不 同 季 节 不 同年 份 的 环境 差异 较 大 ,春季 银川 
平原 土壤 中 的 盐分 不 断 癌 地 表层 汇聚 ,银川 平原 主 
要 为 耕地 ,农民 在 耕种 前 会 对 农田 进行 浸 灌 以 此 来 
降低 农田 土壤 盐 渍 化 程度 ,后 续 研 究 应 该 考虑 不 同 
季节 的 环境 变化 ,选择 多 年 多 季节 数据 进行 预测 模 
型 的 训练 及 验证 。 两 组 输入 参数 对 土壤 盐 渍 化 程 
度 的 预测 精度 不 同 ,在 下 一 步 工作 中 将 对 影响 因子 
与 盐分 指数 进行 组 合 参数 ,进行 区 域 土壤 盐分 的 预 
测 模型 训练 选 优 , 提 高 模型 的 预测 精度 。 男 一 方 
面 , 不 同 地 区 土壤 盐 渍 化 驱动 机 制 不 同 ,本 文 所 建 
立 的 预测 模型 能 否 运 用 到 其 他 地 区 还 需 进行 下 一 
步 的 研究 和 验证 。 


5 结论 


以 银川 平原 为 研究 区 ,将 盐 涡 化 影响 因子 和 盐 
分 指数 作为 输入 参数 构建 盐分 预测 模型 ,对 土壤 0~ 
20 cm 及 20~40 cm 土壤 进行 预测 ,得 到 以 下 结论 : 

(1) 通过 土壤 盐分 预测 模型 对 比 发 现 ,BNN 模 
型 进行 土壤 盐分 预测 效果 最 优 , 影 响 因子 对 0~20 
cm 土壤 盐分 预测 的 影响 较 大 ;盐分 指数 对 20~40 cm 
土壤 盐分 预测 的 影响 较 大 。 对 比分 析 所 选 算法 的 
不 同 变量 建 模 和 验证 效果 ,表明 BNN 模型 是 效果 最 
好 的 预测 模型 ,神经 网 络 的 引入 对 模型 的 训练 有 一 
定 的 优势 。 

(2) 根据 预测 结果 可 知 ,银川 平原 0~20 cm E 
壤 主 要 是 以 非 盐 渍 化 和 轻 度 盐 渍 化 为 主 ,面积 占 比 
分 别 为 33.29% . 31.86% , 重度 盐 渍 化 及 盐 土 面积 
LK 11.59%, 20~40 cm 土壤 中 非 盐 涡 化 和 轻 度 盐 
渍 化 土壤 面积 占 比 分 别 为 50.54% 和 30.07% ,重度 盐 
渍 化 及 盐 土 面 积 占 比 为 7.04%。20~40 cm 土壤 盐 渍 
化 程度 较 0~20 cm 土壤 盐 渍 化 轻 。 
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Abstract: Soil salinization can hinder agricultural development. In this study, the degree of regional soil 
salinization was obtained to provide a theoretical reference for improving agricultural land quality. Using 
Yinchuan Plain of China as the study area with a grid size of 5 kmx5 km, the soil salinity data of 166 sampling 
points at different depths were obtained. Combined with the Landsat 8 OLI image corresponding to the sampling 
time, the salt influence factor and salt index were used as input parameters, respectively, and soil salinity at field 
sampling points was used as output layer parameters. Support vector machine, back propagation neural network, 
and Bayesian neural network (BNN) were established as soil salinity inversion models. The determination 
coefficient and root mean square error of the different models were compared to screen the best model. Finally, 
soil salinization inversion at different depths was performed in the study area. The following results were 
obtained: (1) In the 0-20 cm soil salinity inversion model, the BNN model based on the influence factor variable 
group of salinization was the best, with a coefficient of determination (R°) and root mean square error (RMSE) of 
0.618 and 2.986, respectively; the best inversion result of 20-40 cm soil salinity was the BNN model based on 
the salt index variable group (R°=0.651; RMSE=1.947); the comparative analysis of the modeling and verification 
effects of different variables of the selected algorithms revealed that the BNN model was the best inversion model 
with a better fitting degree than the other two models, and the introduction of a neural network had certain 
advantages in the model construction. (2) Non-salinized and mildly salinized soils were the main soil types in 
Yinchuan Plain. Soil salinization showed a low trend in the south and a high trend in the north. The 20-40 cm 
soil salinization was found to be lighter than the 0-20 cm soil salinization. 
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